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^ ■ Abstract 
0* 

We consider the recently discovered, one parameter family of exactly solvable 
shape invariant potentials which are isospectral to the generalized Poschl- Teller 
potential. By explicitly considering the asymptotic behaviour of the X m Jacobi 
polynomials associated with this system (m = 1,2,3, ...), the scattering amplitude 
\& ' for the one parameter family of potentials is calculated explicitly. 

co ' 

^ '. 1 Introduction 

The ideas of Supersymmetric quantum mechanics (SQM) and shape invariant potential 
(SIP) have not only enriched our understanding of the exactly solvable potentials but have 
helped in substantially increasing the list of exactly solvable potentials [1]. In particular, 
the search for the exactly solvable potentials has been boosted greatly due to the recent 
discovery of exceptional orthogonal polynomials (EOP) (also known as X m Laguerre and 
X m Jacobi polynomials) [4, 2, 3]. Unlike the usual orthogonal polynomials, these EOPs 
start with degree m > 1 and still form a complete orthonormal set with respect to a 
positive definite innerproduct defined over a compact interval. This remarkable work lead 
Quesne [5] to the discovery of two new SIPs (with translation) whose solution is in terms of 
X\ Laguerre and X\ Jacobi polynomials. Soon afterwords, a third SIP (with translation) 
was discovered whose solution is also in terms of Xi Jacobi polynomials [6] . Subsequently, 
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Odake and Sasaki constructed three one parameter family of shape invariant potentials 
(with translation) whose bound state eigenfunctions are in terms of X m Laguerre and X m 
Jacobi polynomials [9]. It is worth reminding here that all of these are isospectral to the 
well known SIPs. 

It is not usually appreciated that unlike the usual SIPs, the newly discovered SIPs 
are explicitly ft dependent. Further, while two out of the three newly discovered SIPs 
have pure bound state spectrum, the third SIP which is isospectral to the generalized 
Poschl- Teller (GPT) potential, has both discrete and continuum spectrum. Recently, we 
have calculated the scattering amplitude for the SIP which is isospectral to GPT poten- 
tial and whose bound state eigenfunction is in terms of X\ Jacobi polynomial [8]. The 
purpose of the present paper is to extend this work to a class of isospectral potentials. In 
particular, in this paper we consider one parameter family of SIPs which are isospectral 
to GPT and whose bound state eigenfunction is given in terms of X m Jacobi polynomials 
(m = 1,2,3...) and obtain the scattering amplitude for this family by considering the 
asymptotic behaviour of the X m exceptional Jacobi polynomials(EOP). 

This paper is organized as follow: In Sec. 2, to set the notation, we briefly review the 
work of Odake and Sasaki [9] regarding the bound state eigenvalues and eigenfunctions 
for the one parameter family of SIPs which are isospectral to GPT and whose bound state 
eigenfunctions are in terms of X m Jacobi polynomials. To motivate our calculation for the 
general case, the scattering amplitude for the potential with the bound state eigenfunction 
in terms of X 2 Jacobi polynomial is discussed in Sec. 3. The most general X m case is 
discussed in Sec. 4. We summarize our conclusions in Sec. 5. 



2 Bound State of Infinitely many shape invariant po- 
tentials 

In this section we essentially set the notation by reviewing the work of Odake and Sasaki 
[9] regarding one parameter family SIPs and the corresponding bound states. We mostly 
adopt their notations in this paper. 

For m > 1, the shape invariant prepotential ui(r; A) which is isospectral to the GPT is 
given by 

f w ( \ i x\ , i £m(cosh2r; X + 5) 
w m (r; A) = cdo{r; A + mo) + log — — — ; < r < oo (1) 

where A = (g,h), h > g > 0, 8 — (1,-1), while £ m (cosh 2r; A) is related with the 
Jacobi polynomial as follows: 

£ m (cosh2r; A) = § )( cosh 2r ) . ( 2 ) 
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It may be noted that the prepotential related with GPT corresponds to m = and is 
given by 

Uo(r; A) = glogsinhr — /i log cosh r . (3) 
The general Hamiltonian corresponding to the prepotential u m (r; A) is given by 

n m (X) = An(A)U m (A) =P 2 + V m (r) , (4) 

where 

p=-id r V m (r)=^ m (r;\) 2 + u'; n (r;\), (5) 

A m (X) = d r - u' m {r; A), A rn {\) ] = -d r -J m {r-\) . (6) 

Here prime on u>i(r; A) denotes derivative with respect to r. The prepotential u m (r; A) 
satisfies the shape invariance condition 

A m (X)A m (Xy = A m (X + 5)A m (X + + S 1 (X + mS) (7) 

This further implies 

u'jr\ A) 2 - u' m (r; A) = u' m (r; A + 5) 2 + J'jr- \ + 5)+£ 1 (X + m5), (8) 

in which 5 is a certain shift of the parameter A. The general form of entire set of discrete 
eigenvalues and corresponding eigenfunctions of H m {X) are obtained by solving, 

/ Hm(X)ilJ m ,v(r; A) = £ m ,„(A)^ m>1/ (r; A) (9) 

The discrete eigenvalues are, 

u-l 

£ m A x ) = + m5) = Y / £i(A + k6 + ttmJ) = Av(h - g - 2m - v) (10) 

k=0 

with v — 0, 1, 2..., z/b — m; z/b = ^ 2 ~ • corres P on ding eigenfunctions are written 

as 

^ m ,u(r; A) = m (r; A)P mjJ ,(cosh 2r; A) 

gU>o(r;A+m£) 

= Ucos/ i 2r;A) P --' COSh2r;A '- <"> 

with 

-P m ,i/(cosh2r; A) = a m ^(r; A) (cosh 2r; A + m<5) + b nhU (r\ A)P,_i(cosh 2r; A + m<5) . (12) 
Here the coefficients a mtU (r; A) and b mtU (r; A) are given by [9] 

2u(-g -h + m- l)£ m _i(r; g,h-2) 



i,u(r; A) = £ m (cosh2r; # + 1, /i - 1) + 



(-g - h + 2m - 2)(g - h + 2is + 2m - 1) 
v(-2h + 4m - 3)£ m _ 2 (r; # + 1, /i - 3) 



(2c/ + 2i/ + !)(-# - /i + 2m - 2) 



, (13) 



, _ (-g - + m - l){2g + 2v + 2m - l^-ifr; g, - 2) 
Wlr,AJ " (2g + 2v + i)(g-h + 2v + 2m-l) ' 1 J 

It is worth noting that the polynomials P m>l ,(r; A) are orthogonal with respect to the 
measure m (r; A) 2 , i.e. 

roc 

/ dr<f) m (r; A) 2 P m ^(r; A)P mi9 (r; A) 

«/ 

= hm, v (g, h)5 vm = h v (g + m,g- m) — + 9 + m ~\ 2 ^ — - — ^ 2 , (15) 

where 

MA)= 2,!(/.-,- 2 ,)r(/.-, + |) • (16) 

It is remarkable that even though the potential related with GPT, i.e. Vqpt = u' (r; A) 2 — 
cu (r; A) is very different from the potential V m = 0J m (r; A) 2 — w m (r; A), the bound state 
spectrum(lO) of the two for any integral m is still the same, however the corresponding 
eigenfunctions are different. Replacing 2r by r, and after using (2), (3) in (11), the bound 
state wave function related to the X m Jacobi polynomial, is given by 



(cosh r _i)^+i/2)( coshr + 1 )|(/3+i/2) ^ 
PiT a - 1 '- /3 - 1) (coshr) 



W(r) = N ^— ,; „ . ^ . ^(coshr) (17) 



where a = g + m- ±,P = -h + m-±, N™ = \2S h - g - 2m+ ^h m , n (g, is the 

normalization constant, Pj+^' m ^ (cosh r) is [y + m) th-degree X m Jacobi Polynomial and 
p(~ a ~ 1 ~P~ 1 \coshr) is usual Jacobi polynomial. 

3 Calculation of scattering amplitude for m=2 (X 2 
Jacobi polynomial) 

The relation between the X m Jacobi polynomial and the usual Jacobi polynomial is given 

by [7] 

^(coshr) = Jftg(.«l,r) = I ^(^r) + ^ - " + & - 



(2m - a + (5 - 2){2v + a + (3) 

_ ^ + m-l)Pi--r ) (coshr) 

(a + is-m + l)(2m-a + p-2) v v ; 

(m - a + g - l)(q + ij p(-«,i9) (coshH p(«,i8) (cosh ^ fl8) 
+ (Q + + i )(2z/ + « + / 3)^- 1 (™sh r)P n _ x (cosh r) (18) 
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where 

P^) (coshr) = 



T(a + u+ 1) 
z/!(I> + /3 + / + l) 



E 

9=1 



z/ \ T(a + /3 + z/ + p + l) / coshr — 1 



V 



(19) 

Using this relation, we have recently calculated the scattering amplitude for the m — 1 
{X 1 Jacobi case) [8]. We now extend that discussion to the m = 2 case. For X 2 Jacobi 
case, we set m = 2 in the above expression, to get 

\w + 2) + („ - - 2 )(o - p - 1)} - 

/ (a-/3-l)(a + /3 + 2) - a + l)(a - /3 - 2)' 



P^f (coshr) 



+ 



4 (/3-a + 2)(/3 + a + 2z/) 

u(/3 -a + l)(a + /3) v(P + l) 



x 



(/3 - a + 2)(/3 + a + 2v) (a + v - l)(/3 - a + 2) 
(/3-a + l)(a + i/) 



Pj a ' /3) (coshr) 
[{a + (3- 2)x + {a + /3)]P^f (coshr) (20) 



The usual Jacobi polynomial Pjf-^ (coshr) can be written in terms of Hypergeometric 
function as : 

p(-»(cosh r) = r( ;;° + ;w + <. + * + 1, i + * iq^) . <») 

z/!l (1 + ctj 2 
To get the scattering state wave functions for this system, two modifications of the bound 
state wavefunctions are required [10]: (i) The second solution of the Schrodinger equation 
which diverges asymptotically and hence had been discarded earlier, must be retained, 
(ii) The discrete level v should be replaced by the wavenumber k such that one gets 
asymptotic behavior in terms of e ±lkr as r — > oo. 
Equation(21) can be written by considering the second solution as, 



P^ (coshr) 



+ Q 



T(v + ct + 1) 
1 — coshr 



1 — prmh t 

C x F(v + a + $ + 1, -u, 1 + a; (22) 



y( v+a +p+i) F ( 1/ + /3 + 1> _ I/ _ a>1 _ a . 



1 — cosh r 



We consider the boundary condition, r — > , i.e.( c ° shr ) — y o, ip v (r) — > finite, the allowed 
solution is 

P^coahr) = r( r + a + ; ) C 1 P(, + a + + 1, -u, 1 + a; iZL^I) (23 ) 
z/!l (1 + a) 2 

where Ci is a constant. In order to compare our results with the previous results [8], we 
use a = B — A — ± , (3 = —B — A — \. Now replacing v by A + ik, we get 



^> <,,,„,,,, - r(S + <Jb + l/2) , F( _ A+ ^_ A _^ jB _ A+1/2; l- coshr. 



P ( ^ (coshr) = 6^ 



(A + iife)!r(S-i4 + l/2)" 



(24) 
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fcoshr) - C. nB + ik- 1/2) 

xF(-A + ik - 1, -A - ik + 1, B - A + 1/2; 1 ~ coshr ) . (25) 

Using Eqs. (24) and (25) in (20) we get Pfcfy (cosh r) = P ( ( ^ +2) (coshr). 
The scattering state wavefunction thus is given by 

ipk(r) = 

n2 (coshr - l)^ B - A )(coshr + lj-^+^P^f/ 1 ^ (coshr) 

fc (2(5 - 1)(2B - 1) cosh 2 r + 2(2P + 1)(2A + 1) coshr + 4A 2 + AA + 25 - l) 

Using the properties of hypergeometric function [11], 

F(a,p, r ,z) = I 1 ~ - r ° r(?)r(7 - «) F( °' 7 -ft"-' 9 + ^ rb> 



and taking the limit r — >■ 00, finally we get the asymptotic form of (26), as 
lim^ fc (r) = 7V 2 - 



2 dP2 A+1 4 ifc 



16 
where 



\PJ 



(28) 



„ (B + ik-3/2){2ik- l)(ab) + 2{B + ik- 1/2) {ed) T(B + ik + 1/2) 

P = tt; : r~r~, — : ; ; d — 



(B + ik-3/2)(2ik-l) ' " (A + ik)\T(B - A + l/2) 3 



T(B-A + l/2)T(-2ik) T(B - A + l/2)T(2ik) 

= —, ; : 7~T ; C 



F(-A - ik)T{B -ik + 1/2) ' F(-A + ik)Y{B + ik + 1/2) ' 



T{B + ik - 1/2) F(B- A+ l/2)T{-2ik + 2) 

d = — ; ^ .„ , „ ; —r—; t = 



(A + ik- 1)\T(B -A + 1/2) ' F(-A -ik+ 1)T(B - ik + 3/2) ' 
The asymptotic behavior for the radial wavefunction (for 1=0) is given by [1] 

lim^ fc (r)~^-[5, = oe ifcr -e- < * r ] (29) 

r^oo 2/c 



From (28) and (29) we get 

v i p , 



S l=0 = ( 2- 4ifc (30) 
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Using the values of P, a, b, c, d and e, we obtain the scattering amplitude for the X 2 
Jacobi case 



Si 



GPT 



{B 2 - (ik - 1/2) 2 } -(B-ik + 1/2) 



{B 2 - 
T(2ik)T(-A 



(ik + 1/2) 2 } - 
-ik)Y(B-ik 



(B + ik + 1/2) 
I- l/2)2~ 4ik 



Y(-A + ik)T(-2ik)T(B + ik + 1/2) 
{B 2 - (ik - 1/2) 2 } -(B-ik + 1/2) 
{B 2 - (ik + 1/2) 2 } - (B + ik + 1/2) _ 



x 



(31) 



4 Calculation of scattering amplitude for X m Jacobi 

case 

We now proceed to generalize this calculation to the X m case (m — 1, 2, 3, ...). Following 
the calculation done above for the X 2 case, using Eqs. (19) and (23) in Eq. (18) and replac- 



(l+m)( coshr ) 



P ( A£Um)( coshr )- Now using P ( 



( ( ^ +m) (coshr) 



ing v — > A + ik we get 

and then taking the limit r — )■ 0, we obtain the asymptotic form of the wave function (17) 
for the Xm case 



-A-2m+ik 



lim ijj k (r) 

r— >oo 

where 



dY(-2B + m- l)r(-2S + 2m - 1)4 

r(-3B + A + 2m - l/2)Y(-2B + m - 1) 



Aik^ikr _|_ ^—ikr 



Q = (ab) 



m(m - 2B - 1)(B + ik - 1/2) 



[B + ik - im + l/2)(2ik - 1) 
while a, b, c, d, and e, are same as in the X 2 case. From (32) and (29) we get 



ac 



-Aik 



(32) 
(33) 

(34) 



Using Q, a, b, c,d and e as given above, we finally have the expression for the scattering 
amplitude in the X m case 



'{B 2 - (ik- l/2) 2 } + (B-ik + 1/2)(1 -m) 
{B 2 - (ik + 1/2) 2 } + (B + ik + 1/2)(1 - m) 
T(2ik)Y(-A - ik)Y(B - ik + l/2)2~ Aik 



-<GPT 
'i=0 



X 



Y(-A + ik)Y(-2ik)Y(B + ik + 1/2) 
{ff 2 - (jfc - 1/2) 2 } + (B - ik + 1/2)(1 - m) 
{5 2 - (ik + 1/2) 2 } + (5 + ik + 1/2)(1 - m) 



(35) 



As expected, in the special case of m — 1 and 2 we get back the expressions for the 
scattering amplitude as obtained in [8] and in Sec. Ill above, thereby providing a powerful 
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check on the calculations. Remarkably, in the limit m — 0, the scattering amplitude as 
given by Eq. (35) reduces to Sj^ 1 , providing a further check on the calculations. It is 
amusing to note that the as one goes from m — 1 to arbitrary integer value, there is 
simply a change by a factor of (1 — m) in the second term in both the numerator and the 
denominator. 

5 Summary 

In this paper we have calculated the scattering amplitude for one parameter family of 
potentials (isospectral to GPT), whose bound state eigenfunctions are given in terms 
of X rn Jacobi polynomials. The bound state eigenvalues and eigenfunctions for these 
potentials were known before [9] . Thus, with the calculation in this paper, one now has a 
complete knowledge about both the bound state spectrum and the scattering amplitude 
for the one parameter family of SIPs which are isospectral to GPT. 
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